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ABSTRACT 

A catalog of common-proper-motion (CPM) companions to stars within 67 pc of the Sun is constructed 
based on the SUPERBLINK proper-motion survey. It contains 1392 CPM pairs with angular separations 
30" < p < 1800", relative proper motion between the two components less than 25 mas yr"', magnitudes 
and colors of the secondaries consistent with those of dwarfs in the {My.V -J) diagram. In addition, we list 21 
candidate white-dwarf CPM companions with separations under 300", about half of which should be physical. 
We estimate a 0.31 fraction of pairs with red-dwarf companions to be physical systems (about 425 objects), 
while the rest (mostly wide pairs) are chance alignments. For each candidate companion, the probability of 
a physical association is evaluated. The distribution of projected separations s of the physical pairs between 
2kAU and 64kAU follows f{s) oc which decreases faster than Opik's law. We find that Solar-mass dwarfs 
have no less than 4.4 ±0.3% companions with separations larger than 2 kAU, or 3.8 ±0.3% per decade of orbital 
separation in the 2 to 16 kAU range. The distribution of mass ratio of those wide companions is approximately 
uniform in the 0.1 <q< 1.0 range, although we observe a dip at <7 2± 0.5 which, if confirmed, could be evidence 
of bimodal distribution of companion masses. New physical CPM companions to two exoplanet host stars are 
discovered. 

Subject headings: stars: binaries 



1. INTRODUCTION 

Very wide binaries play a special role in astro- 
physics. Low binding energies make them sensitive 
probes of Galactic dynamical environment. The very ex- 
istence of wide binaries can constrain properties of the 
dark matter (e. g. Weinberg. Shapiro & WassermannI 119871: 
lYoo. Chaname. & Gould 2004). The typical size of proto- 
stellar cores (~0.1 pc) is comparable to the separation of the 
widest binaries that survive against gravitational perturbations 
in the Galaxy. Most stars form in clusters, in which wide (and 
even relatively closer) binaries are susceptible to disruption; 
the surviving wide pairs provide constraints on the dynami- 
cal e nvironments and p rocesses at the epoch of star forma- 
tion (iParker et al.ll2009h . At larger scales of '-^Ipc we also 
encounter co -moving remnants of nearby clusters and asso- 
ciations (e.g. lZuckerman et al.ll200ll|200 4'). These groups of 
stars, related by common origin, cannot be classified as phys- 
ical systems (binaries) because as far as we know, they do not 
move on stable orbits around each other The transition be- 
tween binaries and co- moving stars is still an unsettled issue. 
IShava&OllingldMl found pairs with projected separations 
of few pc in the Hipparc os catalog, many of which belong 
to known moving groups. ICaballerol(l2010h argues that one 
multiple system with 1-pc separation (member of the Castor 
moving group) could be "hard" enough to remain bound. 

Yet another use of wide binaries is in the calibration of lu- 
minosities and abundances in the low-mass regime, referenc- 
ing to their hotter well-studied primary components. Paral- 
laxes of the bright primaries from, e.g., the Hipparcos cat- 
alog can be assigned to the s econdaries and provide a ccu- 
rate luminosity measurements (iGould & Chanam^l2004l) . A 
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well-characterized primary is particularly useful to calibrate 
the fundamental parameters of M dwarfs, whose complicated 
spectra dominated by molecular absorption bands still resist 
detailed atmospheric modeling. The metallicity scale of M 
dwarfs remains largely dependent on metallicity measure- 
ments o£jroK_priinariesin_TO proper motion (CPM) 
pairs (iRojas-Ayala et al.ll2012h . It is thus important to as- 
semble catalogs of common proper motion pairs where the 
primaries are well-characterised Hipparcos FGK dwarfs, and 
the secondaries are M dwarfs spanning a broad range of colors 
and luminosities. 

The main motivation of this study is to provide a reason- 
ably complete census of very wide companions to nearby 
Solar-type stars. Combined with information on closer 
pairs also available for this sample, it will lead to the un- 
biased multiplicity statisti cs much nee ded as a benchmark 
for star-formation theo ries (iBate Il2009h . The recent study by 
iRaghavan et al.l ( 1201 O f) furnished such statistics for 454 stars 
within 25 pc from the Sun, but a much larger sample is needed 
to study triple and higher-order hierarchies because they com- 
prise only ~10% of all objects. Most solar-type exo-planet 
hosts are also beyond the 25 pc horizon of previous multiplic- 
ity surveys. It is important to search for and identify wide 
companions to known exoplanet hosts, and companions to 
other nearby stars which are still being observed in search of 
exo-planets, i n order to charact erize planet populations in bi- 
nary systems (iRoellet al.ll2012h . 

Primary targets for our search for CPM companions are 
all stars with parall ax tthip > 15 mas selected from the 
Hipparcos-2 catalog (Ivan LeeuwenI l2007l hereafter HIP2). 
The Hipparcos catalog is complete to 67 pc for dwarfs of mid- 
G spectral type or brighter, thus providing a clean distance- 
limited sample for multiplicity statistics. At the same time, 
these nearby stars move fast enough to discriminate CPM 
companions against background stars and are close enough 
to detect intrinsically faint companions. Little is to be gained 
by extending the distance limit further. 

Existing data on CPM binaries is in large part based 
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on the Luyten vis ual searches f or CPM binaries with fi > 
0.1" yr"' (LDS, iLuytenI IT987h . A recent study has re- 
peated the visual search, by blinking t he Digitized Sky Sur- 
vey fields around exo-p lanet host stars dRaghavan et al.ll2006l) 
and stars within 25 pc dRaghavan et al'l 1201 Ol) . A more sys- 
tematic sea rch for companions t o Hipp arcos stars was con- 
ducted by Lepine & Bongiorno| (l2007l hereafter LB2007), 
who used the LSPM-nort h catalog of stars w ith proper 
motions /i > 0.15" yr"' dLepine & Sharal 120051) : the pa- 
per also gives a good overview of prior work. Detecting 
CPM binaries in large da t abases was also demonstrated by 
aname & Gould (2005). Makarov. Z acharias. & Hennessyl 
(|2008), and Dhital et al. (2010), for systems with fainter pri- 
maries. 

In this paper, we expand on the search performed by 
LB2007, and use the full SUPERBLINK all-sky proper mo- 
tion catalog, which is an extension of the published LSPM- 
north catalog. The method for selecting CPM companion can- 
didates is presented in Section |2] Statistics of all companions 
are discussed in Section ID Then we concentrate in Section |5] 
on the statistics of wide physical binaries, namely the distri- 
butions of projected separations and mass ratios. Our conclu- 
sions and discussion of results are presented in Section|6] 

2. IDENTIFICATION OF COMMON PROPER MOTION COMPANIONS 
2.1. SUPERBLINK selection 

The SUPERBLINK survey is based on an automated search 
for moving sources from multi-epoch images of the Digi- 
tized Sky Survey (DSS), notably scans from the first and sec- 
ond epoch Palomar Sky Surveys (POSS-I, POSS-II). Com- 
pared to LSPM-north, the detection threshold in the full SU- 
PERBLINK catalog is significantly lower, /i > 40 mas yr"' 
for stars north of declination -20°, and extends the survey to 
a proper motion limit /x > 150 mas yr"' for all declinations 
south of -20°. The survey is still in pro gress; more infor- 
mation can be found in the recent paper by lLepine & GaidosI 
d2011h . SUPERBLINK is at least 90% complete down to vi- 
sual magnitude V = 19, and has a very low level of false detec- 
tions (< 0.1%) thanks to stringent quality controls including 
visual confirmation of all targets. The V magnitudes of faint 
stars are estimated from the DSS photographic photometry re- 
ported in 'Monet et al.' (^03ft. using the procedure described 
in Lepine & Shara (2005). 

We select as possible primaries of CPM pairs all stars from 
HIP2 with parallax greater than 15 mas and proper motion 
larger than the SUPERBLINK limit (/x > 40 mas yr"' for stars 
north of declination -20° and /i > 150 mas yr" otherwise). 
Potential CPM companions to these stars were selected from 
SUPERBLINK by the following three criteria. 

1 . Separation of the companion p should be in the interval 
30" < p < 1800". Closer companions are affected by 
the oreols of primary stars in DSS photographic images, 
introducing selection against faint companions. Most 
pairs wider than 300" (projected separation 20kAU at 
67 pc) are chance alignments (optical), but including 
them in the catalog helps us to evaluate the proportion 
of true (physical) wide binaries. 

2. Proper motion difference with the primary component 
Ap, should be less than 25 mas yr"'. Typical proper 
motion errors in SUPERBLINK are 8 mas yr"', so this 
cutoff allows for larger errors or some additional proper 
motion difference caused by motion in hierarchical sub- 
systems. 
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Fig. 1 . — Color magnitude diagram of the common proper motion secon- 
daries, assuming they have the same parallax as the (Hipparcos) primaries. 
Physical companions are expected to fall either close to the main sequence 
(red) or the white dwarf sequence (blue); a scatter of ±1.5 mag is adopted 
for the selection (shown). Stars falling outside the bounds are assumed to be 
chance alignments of unrelated stars. Top: diagram for the 608 pairs with 
separations 30" < p < 300". Bottom: diagram for the 240 pairs with separa- 
tions 1470" < p < 1500". which cover a similar search area. Pairs in the first 
sample are more numerous and tend to concentrate within the expected main 
sequence and white dwarf sequence, suggesting that the f?i400 extra pairs in 
the first distribution are physically related. 



3. Color and magnitude of the secondary should be con- 
sistent with a main sequence star or a white dwarf at the 
same distance as the primary. This criterion effectively 
requires that the photometric distance of the secondary 
as estimated from {V,V -J) be consistent with the par- 
allax distance of the alleged Hipparcos primary. 

The first two criteria identify 130882 possible pairings, the 
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vast majority of which are simple chance alignments. For the 
third criterion, we first calculate the absolute magnitudes of 
the CPM companions {[Mv]5ec) assuming that the secondary 
is at the same distance as the Hipparcos primary: 



200 



[My ] SEC = VsEC + 5 log(7rpRi ) + 5 , 



(1) 



where ttpri is the Hipparcos parallax of the primary. We then 
use the {My ,V-J ) color magnitude relationship as calibrated 
in (iLepind 120051) for main sequence (MS) stars in the SU- 
PERBLINK survey. These define photometric absolute mag- 
nitudes [Mv]ms estimated from V-/: 

[Mv]Ms=0.08 + 3.89(y-7) [y-7<1.5] 

[MvIms =2.78 + 2.09(y-y) [1.5<V-J< 3.0] 

[Mv]ms = 1.49 + 2.52(y-7) [3.0 <V-J< 4.0] 

[Mv]ms =2.17 + 2.35(y-y) [4.0 <V-J< 5.0] 

[Mv]MS=4.47+1.89(y-7) [5.0 < y-7 < 9.0]. (2) 

In addition, we introduce the following simple color- 
magnitude relationship which we use to define absolute mag- 
nitudes for white dwarfs (WD) based on y-7: 

[Mv]wD= 12.5 + 3.75(y-7) [-i<y-y<0] 

[Mv]wD = 12.5+1. 89(y-y) [0<y-7<3]. (3) 

We then select only pairs whose [Mv]5ec agree to within 
1.5 magnitudes of either photometrically estimated abso- 
lute magnitudes [Mv]ms or [Mv]wd- For 124 companions 
without 2MASS photometry the V - J color is estimated 
from the (fe,r,i) phot ographic magnitudes as described by 
iLepine & Sharal 120051) . 

The color-magnitude selection is illustrated in Figure [T] 
The top panel shows [Mv]sec as a function of y-7 for all 
608 pairs with angular separations 30" < p < 300", with the 
selection range for MS stars in red, and the selection range for 
WDs shown in blue. The bottom panel in Figure [T] shows the 
same distribution but for 240 pairs with angular separations in 
the 1470" < p < 1500" range. Note that the two ranges cover 
the same area on the sky, so if all CPM pairs were the result 
of chance alignments, one would expect the two samples to 
have similar numbers. Instead we find an excess of 368 pairs 
at closer separations. In addition, the color-magnitude distri- 
bution is very different in the two samples. Most pairs with 
1470" <p< 1500" fall outside of the MS or WD bands, while 
the majority of the close pairs have colors consistent with true 
companions. This demonstrates that ^ 400 CPM pairs in our 
catalog are physical, binary/multiple systems. 

The color-magnitude selection reduces the catalog to only 
2 078 pairs with consistent colors. A histogram of the number 
of pairs as a function of their angular separation p is shown 
in Figure ID The distribution for pairs with consistent paral- 
lax/photometric distances (bottom panel) is compared to the 
distribution of pairs where the primary and secondary do not 
have colors consistent with physical association (top panel). 
The latter distribution shows the trend expected from chance 
alignments, i.e. a linear increase with p proportional to the 
increase in the area of sky; this is modeled by the red line. 
The distribution of stars with consistent MS or WD color- 
magnitude, however, shows a marked excess of pairs with 
p < 200", which again points to the existence of physical 
pairs. At larger angular separations p > 450", the distribu- 
tion shows the linear increase expected from chance align- 
ments; the red line shows the predicted trend based on the 
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Fig. 2. — Histogram of the number of pairs found as a function of angular 
separation p. Top: distribution for pairs assumed to be chance alignments, 
based on the inconsistency in the parallax of the primary and photometry 
of the secondary, i.e. stars falling outside of the MS or WD sequences (see 
Fig. [T). Bottom: distribution for pairs in which the secondary would have 
(V,V—J) consistent with a main sequence star or white dwarf. In either case, 
the red line shows the distribution expected from chance alignments (linear 
increase with p). Pairs in which the color/magnitude of the secondary is 
consistent with a physical companion show a net excess at p < 200" . 



number of pairs with p > 1000". One notices that while a 
significant number of physical pairs are undoubtedly detected 
at p < 450", this sub-sample still suffers from contamination 
from chance alignments. 

We note that some companions share common proper mo- 
tion with more than one Hipparcos star. Likewise some Hip- 
parcos stars have several CPM companions. Some are gen- 
uine multiple systems (triples), but in most cases they are 
more likely to be the result of chance alignments, e.g. the 
alignment of a foreground/background proper motion star 
with a true physical pair. In any case, all possible pairings 
are retained in our list. 

2.2. Treatment of the White Dwarf CPM Companions 

Magnitudes and colors of 386 companions place them in 
the WD band in Figure [T] They are found in smaller num- 
bers than MS candidates, and are significantly more affected 
by contamination. Indeed, only 21 pairs have p < 300", and 
some of those are expected to be chance alignments. We eval- 
uated the number of physical WD companions within some 
separation p by their excess over the linear distribution, as in 
Figure |2l and found that the number of likely physical pairs 
with WD companions does not increase beyond p > 200", 
reaching an excess of around 10 pairs. 

Since the vast majority of the stars falling in the WD region 
are thus chance alignments, we removed form the catalog all 
companions in the WD band with p > 300", keeping only the 
21 most likely pairs. 

2.3. Exclusion of nearby cluster stars 

The selection of CPM pairs is sensitive to faint mem- 
bers of the Hyades and Pleiades clusters, whose members 
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have proper motions within the /i > 40 mas yr"' limit of 
the SUPERBLINK catalog, and are thus generally detected. 
Some nearby young associations and m oving groups also hav e 
members detected in SUPERBLINK jSchlieder et al.ll2012l) . 
Faint cluster stars share a proper motion with the brighter 
members, which are also Hipparcos sources. Cluster stars 
have photometric distances consistent with the parallax of the 
brighter members, so they will generally meet our selection 
criteria for CPM pairs with colors and magnitudes consistent 
with physical companions. 

We wish, however, to exclude those pairings because they 
are not physical binaries. We thus eliminate 311 pairs where 
the primary stars are known members of clusters and asso- 
ciatio ns according to the XHIP catalog (I Anderson & Francis I 
120121) . Of those, 230 belong to the Hyades. This leaves 1392 
CPM companions with MS colors and 21 companions with 
WD colors in the catalog. 

2.4. Multiple stars 

In the catalog compilation and analysis presented below we 
assume implicitly that both the primary and wide secondary 
companion are themselves single stars, i.e. there are no un- 
resolved companions. In reality a substantial but still poorly 
known fraction of the com ponents in wide systems are them- 
selve s close binaries (e.g. iMakarov. Zacharias. & Hennessvl 
|2008|) . Internal orbital motion can contribute additional com- 
ponents to the proper motions and increase A/i with the pri- 
mary. Analysis of a strometric binaries a mong FG dwarfs 
within 67 pc listed by iMakarov & KaplanI (l2005h shows that 
for 86% of them the difference between long-term and Hip- 
parcos proper motion is less than 25 mas yr"'. Therefore, our 
catalog may be slightly biased against triple systems. 

To verify this, we simulated a realistic population of sub- 
systems around primary targets with periods from 1 yr to 
lOOOyr susceptible to create A/i by photo-center motion dur- 
ing the 3.2-yr long Hipparcos mission. Maximum effect is 
produced at periods from 4 yr to 40 yr, but the mean and me- 
dian A/i do not exceed 6 mas yr"' at 50 pc distance; in 90% 
of cases A/t < 12 mas yr"' for all periods. According to these 
simulations, the fraction of multiple systems rejected by the 
A/i < 25 mas yr"' threshold should be sign ificantly smaller 
than the 14% implied by the paper of Makar ov & KaplanI 
((2005). The discrepancy is alleviated by including in the sim- 
ulation some ma ssive companions that are close binaries and 
cause larger A/i dTokovinin et al.ll2012l) . 

Multiple systems can also be missed because of inaccu- 
rate or distorted photometry. Visual magnitudes in the SU- 
PERBLINK catalog are largely based on photographic mea- 
surements, and sometimes have large errors, especially if the 
star is blended with another source, such as in the vicinity of 
a very bright primary or in a dense field. The composite mag- 
nitudes of unresolved sub-systems also bias photometric dis- 
tance estimates which might yield to rejection based on color 
and magnitude - although the 1 .5 magnitude selection range 
(see Section 2.1) mitigates this effect. In any case, this may 
be another non-negligible, but poorly quantified, incomplete- 
ness factor for hierarchical systems with CPM companions. 

3. CATALOG DESCRIPTION 

The complete catalog of 1413 CPM pairs is listed in Ta- 
ble 1, the full version of which is available electronically. 
The presentation has some similarity with LB2007. The 
first columns contain information on the primary target: (1) 
Hipparcos number HIPpRi, (2) SUPERBLINK identification. 



(3,4) equatorial ICRS coordinates at the 2000.0 epoch in de- 
grees, (5,6) proper motions in RA and Dec in arcsec yr"', (7) 
V magnitudes derived from the DSS for faint stars or taken 
from Tycho-2 for bright ones, (8-11) V -J color and J,H,K 
magnitudes from t he 2MASS All-Sky Point Source Catalog 
(ICutriet al. Il2003h . Then follow the parallax (12), its error 
(13), and the parallax code (14), namely PLXl for trigono- 
metric parallax from HIP2. The following columns list (15) 
separation p in arcseconds, (16) proper motion difference with 
the primary target A/i in mas yr"', and (17) estimated prob- 
ability of physical association Pphys (Section [4. 4l i. Compan- 
ions in the WD band (see Section 2.2) are distinguished by 
^phys = -1 and listed at the end. 

The remaining columns of Table 1 contain data on the 
CPM companions. Columns (18) to (31) give same infor- 
mation for companion as columns (1) to (14) for the pri- 
mary. The HIPsEC number in (15) is 999999 if not available. 
Similarly, missing magnitudes are listed as 99. Photomet- 
ric parallaxes in column (29) estimated from the My ,V -J 
relation have code PHOT in column (31), their errors are 
set to zero. Parallax of some companions is taken from 
the hte rature and marked by the following codes: PLX2 
dyanAltena et al. 1995)'', PLX3 from the NSTARS database^ 
and PLX4 (Mvers et al. 2002)^. Trigonometric parallaxes of 
two WDs ( McCook& Sion .,19991) have code PLX5, the pho- 
tometric parallaxes of WDs are calculated from (O. 

We emphasize again that this catalog contains all CPM pairs 
identified using the selection method described in §2, and in- 
cludes not only physical binaries but also a significant number 
of chance alignments. One should therefore not assume that 
the catalog lists true companions of the Hipparcos stars, and 
the user should keep an eye on column 17, which provides a 
probability for the pair to be physical. We describe below our 
statistical method to assign probabilities of physical associa- 
tion for all the CPM pairs. 

The 21 candidate WD companions are highlighted in Ta- 
ble |2] extracted from the main catalog. The Hipparcos num- 
ber of the primary component HIPpRi is listed in column (1), 
followed by the equatorial coordinates of the companion for 
2000.0 in hexadecimal format in columns (2) and (3). Subse- 
quent columns contain separation (4), A/i (5), proper motion 
of the primary component (6), its parallax (7), Vsec magni- 
tude of the companion (8), and its V -J color or its estimate 
(9). The last column (10) gives the identifications of 5 previ- 
ously known WDs found in SIMBAD. Another 6 promising 
WD candidates are marked as wd? in this column on the basis 
of their large proper motion or small separation. The compan- 
ion to HIP 1 18010 is also Usted in LB2007 . HIP 60081 hosts 
a planet with 48 d period (ISchneidedl201 lb . 

4. STATISTICAL PROPERTIES OF THE CATALOG 

Half of our targets have /i < 87 mas yr"' and 3/4 have /i < 
150 mas yr"'. Therefore the extension of the SUPERBLINK 
proper motion limit to 40 mas yr"' is critical for this work. On 
the other hand, the catalog is still incomplete atS < -20°. The 
distribution of companions on the sky is uniform, but with 
fewer stars south of S = -20° due to the larger proper mo- 
tion cutoff. Interestingly, there is no obvious concentration 
towards the Galactic plane. Although the density of back- 
ground stars strongly depends on the Galactic latitude, they 

^ http : //vizi er. cfa. h arvard. edu/viz-bin/Vizie R?-sou 
^ http : / / nst ars ■ nau ■ ed u/ nau_nstars/ index . htm 
http: //cdsarc.u-strasbg.fr/viz-bin/CatTV/lOgl 
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Fig. 3. — Proper motion difference Afi vs. angular separation p for the 
1 392 common proper motion pairs with MS secondaries in our catalog. Pairs 
with short separations are dominated by physical systems, whereas pairs with 
large separations are mainly chance alignments. Vertical dashed Hnes mark 
separations of 200" and 1200". 
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Fig. 4. — Cumulative distributions of Afi for close (p < 200", N = 315, 
dashed line) and wide (p > 1200", N = 552, full line) CPM pairs. The crosses 
show analytical fits. Close pairs show a distribution which decreases in Ap, 
consistent with physical systems, while wider pairs show one that is increas- 
ing in Afi, consistent with chance aUgimients. 



mostly have sr nall u and thus are not in cluded as companion 
candidates (see lLepine & Gaidosll201 ih . 

Selection of primary targets with relatively large proper mo- 
tion naturally favors nearby stars. Median parallax of the cat- 
alog entries is 18.1 mas. The number of FG dwarfs within 
distance d is proportional to li^ ' in our catalog, whereas it 
is close to in Hipparcos; the difference is caused by the 
proper motion selection. In this Section we do not consider 
the 21 WD candidate companions and analyze only the 1392 
MS companions. 

The parallax distributions for close and wide companions 
in our catalog are slightly different; the median parallax is 
22.1 mas and 17.1 mas for p < 300" and p > 1200", respec- 
tively. Distant stars have, on average, a smaller proper motion 
and a larger number of optical companions, producing this 
correlation between tthip and p. 



4.1. Proper motion difference 
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Fig. 5. — Visual magnitude of the secondaries VsEC as a function of the 
angular separation p of the CPM pair. The secondaries in close (p < 100") 
parrs, which are largely physical systems, have a relatively uniform magni- 
tude distribution. Secondaries in wider (and mostly change alignment) pairs 
show a magnitude gap at VsEC ~ 13 which can be explained by a combination 
of selection effects. 



Figure [3] plots proper motion difference A/i for the pairs 
with main sequence secondaries against their angular sepa- 
ration p. The distribution reveals statistical differences be- 
tween physical (close) and optical (wide) companions. Com- 
panions with p < 200" have uniform distribution of log p and 
mostly A/i < 15 mas yr"'. In contrast, wide companions with 
p > 1200" tend to have larger A/i. The difference of the A/i 
statistics between these two groups is illustrated in Figure |4l 
while 56% of close companions have A/i < lOmas yr"', this 
cutoff contains only 19% of wide companions. Cumulative 
distribution for wide companions can be fitted by a power law: 



Fop,(A/i) « (Ap/25) 



1.78 



(4) 



If the proper motion vectors of background stars were 
distributed uniformly, one would expect a quadratic law 
(LB2007), but the actual power index is slightly less than two; 
this is because nearby stars have systemic motions relative to 
the Sun and their proper motions thus show local correlations 
on the sky. The analytical model for close (physical) compan- 
ions is: 



Fphys(A//)«[(A/i-2)/22]" 



(5) 



These two distributions overlap substantially, preventing us 
from using a smaller cutoff in A/i. Both distributions extend 
beyond the 25 mas yr"' cutoff (the cumulative curves have 
non-zero slope near the cutoff), suggesting that the catalog 
misses some physical companions with A/i > 25 mas yr"', 
e.g. wide triples with motion in inner sub-systems. Physical 
companions with Ap > 25 mas yr"' are indeed identified in 
LB2007. However, the large number of chance alignments 
with A/i > 25 mas yr"' places very low probabilities on any 
pair being physical and would require significant follow-up 
resources for triage and confirmation. For now, we prefer to 
adhere to this (somewhat arbitrary) upper limit in Ap. 

4.2. Magnitudes and colors 

In our catalog of CPM binaries, secondaries in wide pairs 
tend to be fainter than those in closer pairs (Figure|5]l. Detec- 
tion of faint close companions is affected by the oreols around 
bright primary targets (LB2007), but at p > 30" it seems to be 
complete at least down to Vsec = 17"'. 
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Fig. 6. — Cumulative distributions of AA/v for close (p < 200", A' = 315, FiG. 7. — Cumulative distribution of separations (full line) and its model 

dashed line) and wide (p > 1200", A' = 552, full line) companions. Analytical (dashed line). The dotted line corresponds to physical companions, 

models are plotted as crosses. 



There appears to be a deficiency of companions with 
VsEC from 12™ to 14"'. Bright candidates with Vsec < 12 
were selected from the Tycho-2, while fainter candidates 
are identified directly from the DSS with SUPERBLINK. 
There is, however, an overlap of a couple magnitudes be- 
tween the Tycho-2 catalog and SUPERBLINK identifications. 
|L6pine & Shara (2005) argued that the merged Tycho-2 + SU- 
PERBLINK catalog should be 98% complete in this overlap- 
ping range of magnitudes, based on the high rate of recovery 
by SUPERBLINK of the faintest of the Tycho-2 stars. How- 
ever, the Vsec magnitudes of faint companions are still derived 
from the DSS scans of photographic plates, with estimated er- 
rors of ±0.5™ or larger At the bright end, saturation of the 
stars on the plates may lead to systematic errors, and a slight 
over-estimation of V could create the discontinuity of the Vsec 
distribution seen in Figure |5] 

There is, however, an alternative explanation. One can see 
from Figure |5] that the gap is most apparent in the very wide 
secondaries, which are largely dominated with chance align- 
ments, and much less obvious in the close (p < 100") pairs. 
While any field star could potentially be selected as CPM sec- 
ondary based on proper motion alone, our color-magnitude re- 
quirement eliminates most of the foreground and background 
contaminants, because their colors and magnitude would not 
fit the expected MS at the distance of the star. Exceptions 
to this rule are background giant and subgiant stars, which 
generally cannot be distinguished from foreground MS stars 
based on color and magnitude. Effectively, this causes an 
excess of (chance alignment) wide CPM secondaries in the 
7 < Vsec < 12 magnitude range; fainter subgiants are not 
found in the SUPERBLINK catalog because of the distance 
bias introduced by the high proper motion cut. This over- 
selection of bright giants/subgiants conspires with the rising 
luminosity function of field stars to create the magnitude gap 
for wide companions. 

A color-magnitude diagram (CMD) of the companion can- 
didates is plotted in Figure [T] There is no clear separation 
between optical and physical companions in this CMD, but 
the statistics of the deviation from the nominal MS AMy = 
[Mv]sec-[Mv]ms (see Section 2. 1) is different between those 
two groups. Most close (physical) companions have AMy ^ 
0, while wide companions tend to be fainter and mostly have 
positive AMy. Cumulative distributions of AMy are com- 



pared in Figure |6] The cumulative distribution of AMy for 
wide companions is approximated by a cubic polynomial in 
AMv with coefficients (0.2052,0.2342,0.1210,0.0468). The 
AMy distribution of close companions is however normal 
with a dispersion of 0.575'" and zero average. Our cutoff of 
±1.5"' (2.6(7) rejects only 1% of this normal distribution. This 
indicates that our color-magnitude selection misses very few 
MS companions. A reduction of the cutoff would likely re- 
duce the contamination by optical pairs, but would result in a 
lower completeness for physical systems. 

4.3. Distribution of angular separations 

We expect that the number of optical (i.e. chance align- 
ment) companions closer than p should be proportional to p^. 
Physical companions, on the other hand, are exp ected to b e 
distributed almost uniformly in log p (Opik's law, OpikI 19241) . 
This allows us to model the cumulative distribution of separa- 
tions as a sum of physical and optical populations, 

F(p) = ax+bx^ + c(p/p2)\ (6) 

where: 

X = logip/pi)/ l0g(p2/ Pi)- 

Here pi = 30" and p2 = 1800" are the separation hmits of the 
catalog. The variable x maps this interval logarithmically into 
(0, 1), and the parameters a and b describe the distribution of 
physical companions linear in x (instead of uniform distribu- 
tion we allow some slope b), while c = I -a- bis the fraction 
of optical companions. This simple model with only two free 
parameters fits well the observed distribution of separations in 
Figure|2]with a = 0.600 and b = -0.294. The estimated frac- 
tion of optical companions in our catalog is therefore c = 0.69, 
so it appears to contain about 425 physical wide pairs. The 
same model predicts that 92.7% of companions with p < 300" 
are physical. 

4.4. Probability of physical association 

We have modeled the distributions of two astrometric pa- 
rameters, angular separation p and proper motion difference 
A/i, by fitting analytical formulae to the cumulative his- 
tograms (eqs. We also modeled the distribution of the 
deviation from the main sequence AMy- Probability den- 
sity functions for physical (close) and optical (wide) com- 
panions are readily obtained by differentiating these formu- 
lae. We assume here statistical independence of variables 



Wide companions to nearby Hipparcos stars 



7 



p, A/i, and AMy, which is not quite true owing to correla- 
tions between these parameters, but still a reasonable approx- 
imation. Neglecting the correlations, the three-dimensional 
probability density functions for physical and optical com- 
panions are products of their one-dimensional distributions, 
f(p, Ap, AMv) = fp(p)fA^(Ap)fAMv(AMv). 

Each companion can be either physical or optical. So the 
probability of physical association can be estimated as 

PphysCp, A^, AMv) = /phys/(/phys + /opt) (7) 

provided that the distributions /phys and /opt are normalized in 
the same way. We list estimates of fphys in the catalog, but 
do not use them in the statistical analysis presented below. 
The values of Pphys cluster near zero and one, discriminating 
well most candidates. Only 133 companions have interme- 
diate fphys between 0.2 and 0.8, most of them with separa- 
tions of few hundred arcseconds. We checked the estimates 
of Pphys on 175 companions with trigonometric parallaxes of 
±5 mas or better precision and found that for all pairs with 
^phys > 0.5 the parallaxes of primary and secondary compan- 
ions match within errors, thus confirming the physical nature 
of those pairs. However, 39 wide pairs with Pphys < 0.5 (out 
of 70) also have parallaxes that match within 5 mas and are 
probably physical. Our estimates of Pphys based on empirical 
data modeling are therefore useful as a guidance but should 
n ot be taken too literally. 

iLepine & Bongiornol ( l2007h also included the mean proper 
motion of the pair p in their statistical model and showed 
that the probability of optical companions is proportional to 
/i"^ i.e. that CPM pairs with larger proper motions are more 
likely to be physical. This analysis, however, neglected the 
role of color-magnitude selection in eliminating background 
contaminants, here represented by /amv ■ In addition, LB2007 
used A^ oc p~' as a cutoff, which introduced incompleteness 
at large separations by restricting the range of Ap values. For 
these reasons we decided not to include p in the statistical 
criteria used to separate physical and optical companions, in- 
stead relying on /amv to eliminate most optical pairs at large 
angular separations. 

To summarize, the current catalog of CPM companions to 
Hipparcos stars contains a mixture of physical and optical 
companions with different statistical properties such as dis- 
tributions of separations, Ap, magnitudes, and colors. We de- 
velop a model for these distributions which we use to evaluate 
the probability of membership in either of the two groups. Ad- 
ditional parameters, such as total proper motion p, could be 
used to refine those probability estimates, but are not applied 
here for simplicity and uniformity. 

5. STATISTICS OF WIDE BINARIES 

5.1. Masses and mass ratios 

The masses of the primary and secondary components are 
evaluated her e from their absolute magn itudes Mk using the 
relation from iHenrv & McCarthy I ( 11993 ) for dwarfs and as- 
suming that both companions are single (although binarity has 
little effect on estimated mass). The mass-magnitude relation 
in the K band is less sensitive to age and metallicity, compared 
to optical magnitudes. Evolved stars do not obey the relation 
for dwarfs, and as a result are assigned masses above 1.5Mq 
owing to their high luminosity. In any case, 87% of primary 
components have estimated masses between 0.6 and 1.5M0, 
i.e. are nearby solar-type dwarfs. The narrow range of masses 
for the primaries is of course achieved by design, from a com- 
bination of the 67 pc distance limit and the magnitude limit 



1 .4 7- 
1 .2 - 
1 .0 '- 
0.8 '- 
0.6 7 
0.4 7 

0.2 

0.0 
0.0 



q = 1 
q = 0.1 



+\ ++ + 



0.5 



1 .0 



1 .5 
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of the Hipparcos catalog itself. Massive stars are rare in the 
vicinity of the Sun, and most <O.6M0 dwarfs in the 67 pc 
volume are below the Hipparcos magnitude limit. The sec- 
ondary components, in contrast, span a range of masses from 
0. 1 Mq up to the mass of the primary component, consistent 
with our broad detection of secondary stars down to the V = 20 
magnitude limit of the SUPERBLINK survey. 

A comparison of the estimated masses of the primary and 
secondary components is shown in Figure [8] for the 403 pairs 
with Pphys > 0.5. The majority of mass ratios q = Msec/A^pri 
are comprised between 0. 1 and 1 ; we set ^ = 1 in a few cases 
where Msec > A/pRi (which happens because the V magni- 
tudes are used to define which star is the primary/secondary 
whereas the K magnitudes are used to estimate their masses). 

The scarcity of very low-mass companions is obvious in 
the CMD (Figure[Tll by the scarcity of points with V-J > 5.5 
or spectral types later than M6V. At 67 pc (distance modu- 
lus 4.12) a y = 19 star near the SUPERBLINK magnitude 
limit corresponds to a 0. 12Mq d warf with y-7 = 5.2, accord- 
ing to the standard relations of iHenrv & McCarthy I (Il993h . 
Therefore, the catalog is incomplete below Msec ~ O.12M0, 
although CPM companions of lower mass are found around 
some of the nearest primaries in the catalog. 

The mass ratio of wide binaries has approximately uni- 
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form distribution (Figure |9|l. Excluding the first bin af- 
fected by incompleteness, we fit the distribution by a power 
law f(q) PC (7"'^ with /? = 0.07. This matches the result of 
iRaghavan et al.l (12010 ) who found uniformly distributed mass 
ratio for all binaries independently of their period, and the re- 
sult of Tokovinin (201 1) who studied the mass-ratio distribu- 
tion of companions within 20" and also found it to be uniform. 
Hierarchical multiples are neglected by this analysis. 

However, we observe a prominent local minimum around 
q ^ 0.5, which is apparent in Figure |9] and corresponds to 
companions of ^ Q.5Mq mass with [Mv]sec ^ 10 and V-J ^ 
3.4. At a typical distance of 50 pc such companions have 
VsEC ^ 13.5. As noted above, SUPERBLINK is expected 
to be complete in this magnitude range. However, system- 
atic errors in the Vsec magnitudes may affect the position in 
the CMD (Figure 1), l eading to a r ejection of some physical 
companions. Recently R ical (1201 2|) found CPM companions 
to exo-planet hosts with high proper motion; one of those, 
HIP 112414(p = 51", Vsec = 14.1, Ai = 214mas yr"') is over- 
looked in our catalog. Yet, we see only few such potentially 
rejected companions on the right side of the MS band in Fig- 
ure[ri which suggests that such rejections should be relatively 
unfrequent, and cannot explain the observed deficit at ^ '--^ 0.5. 

Companion masses estimated from the K^ec magnitudes 
are not affected by an y systematic errors of Vsec- However, 
the standard relation of iHenrv & McCarthy I (Il993h used here 
changes slope at Mk = 5.9, contributing to the deficit of esti- 
mated masses around ^ O.5M0. Use of an alternative smooth 
relation between mass and Mk decreases the depth of the min- 
imum around q k, 0.5, without however getting rid of it en- 
tirely. 

Should the minimum be real, it might indicate a bimodal 
distribution in the masses of the wide secondaries. A popu- 
lation of roughly equal-mass systems, combined with another 
where the secondaries have a mass distribution similar to field 
stars, could produce a bimodal distribution of mass ratios such 
as the one suggested here. Because of the relatively small 
number of systems in the present study, this apparent deficit 
?X q K, 0.5 should be confirmed before it can be accepted as 
real. 

5.2. Projected separations 

Projected orbital separations s for our wide systems can be 
calculated from s = p/nuip- The separations are statistically 
related to the orbital semi-major axis a. Simulations show that 
in most cases the ratio s /a is comprised between 0.5 and 2 and 
the median of i/a is close to 1 .0 (it depends slightly on the ec- 
centricity distribution). We plot the distribution of orbital sep- 
arations against the distance of the pair from the Sun in Fig- 
ure [TO] The range of projected orbital separations over which 
our survey is sensitive is defined by the lower and upper limits 
of angular separations, which are noted by the dashed lines. 
The majority of the pairs have 33pc < d < 67pc; over that 
range our census is complete for systems with projected or- 
bital separations 2kAU <s < 64kAU. In that regime, the dis- 
tribution of projected angular separations appears to be uni- 
form over log(5) up to a projected separation of lOkAU, be- 
yond which the density increases significantly. This increase 
in density marks the increasing contamination from optical 
pairs (i.e. chance alignments). 

For evaluating the companion frequency, we define the ref- 
erence sample of Hipparcos stars that fulfill the selection cri- 
teria of our catalog on distance (tthip > 15 mas) and proper 
motion (/i >40mas yr"' north of -20° and /i > 150 mas yr"' 
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Fig. 10. — Distribution of projected orbital separations plotted against 
distance from the Sun. The projected separations are calculated assuming 
the secondary is at the same distance from the Sun as the Hipparcos primary. 
Angular separation limits of the survey are indicated by the dashed line. Our 
census is complete for physical pairs with .v > 2 000 AU, but contamination 
from optical pairs becomes significant at .v > 20000 AU. The boxes show the 
range of the five bins in separation and distance used in our estimate of the 
distribution of orbital separations. 



Otherwise). We evaluate their masses by the same method 
as for the catalog targets using .^T-band magnitudes from 
XHIP d Anderson & Francis 1120 121) : for a small subset without 
2MASS photometry, V magnitudes were used instead. The 
reference sample is not complete within its distance limit. It 
only serves to relate the number of CPM companions in our 
catalog to the total number of primary stars with similar char- 
acteristics. 

In the following, we apply additional selection criteria on 
tthip, /i, and mass of primary components to both our catalog 
and reference sample. We also limit the analysis to primary 
stars in the mass range 0.7 Mq <M<1.4M0, which includes 
77% of the pairs in our catalog. 

To derive the distribution of projected separations, we need 
to account for the optical companions which contaminate the 
sample. We cannot use here the estimates of Pphys because 
they depend on p and hence bias the result. Considering a 
restricted range of parallax values, we estimate the number 
of optical companions in five bins of projected orbital sepa- 
rations (see Figure [TOb. by translating their limits (si,S2) into 
angular separation limits (pi , pa) using the average parallax of 



the selected stars. Then A^opt =^opt(/02~P?)' where the param- 
eter Aopt is evaluated from the number of companions between 
1200"and 1800", assuming they are all optical in that range 
(a slight over-estimate). The estimate of Aopt relies mostly on 
companions with s > 64 kAU, outside the range of separations 
studied here. 

The procedure to derive the separation distribution is thus 
as follows. A subset of the catalog satisfying some additional 
criteria is chosen and the number of reference targets A^ref sat- 
isfying the same criteria is found. The histogram of projected 
separations A^, is constructed and the estimated number of op- 
tical pair contaminants is subtracted from each bin / to get 
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histogram, dashed line - estimated optical companions, dotted line and pluses 
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A^i.phys = Ni-Ni,opt- The companion frequency is A^/.phy,s/A^ref- 
We ignore the tiny incompleteness of the fifth (32-64) kAU 
bin. 

The top panel of Figure [TT] shows the result for all stars 
with solar-type primaries (masses from 0.7 Mq to IAMq). 
This selects 518 companions, 227 of which should be phys- 
ical, while A^ref = 5196. The total frequency of companions 
between 2 kAU and 64 kAU is therefore 4.4 ±0.3%. We do 
not find any strong dependence of companion fraction on the 
mass of the primary. The total binary fraction, however, crit- 
ically depends on the subtraction of optical pairs, especially 
in the the last two bins which are significantly affected by the 
contamination. The first three bins contain 180 physical com- 
panions, or 3.8±0.3% per decade of separation, and are not 
much affected by the contamination. 

Our catalog is incomplete for companions of very low 
masses {q < 0.15). The magnitude of this incompleteness 
depends on the actual distribution of q: if the distribution is 
uniform all the way down to q = 0, then the catalog misses 
only 15% of low-mass companions. However, the number 
of sub-stellar companions is probably small (the so-called 
brown-dwarf desert) and the incompleteness at low q is corre- 
spondingly less. Wide triples may also be under-represented 
due to our color-magnitude selection (see 32.4l i. though the 
losses are probably modest. Assuming that our census of wide 
companions is 80% complete, their true frequency would be 
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5.5%. This value corroborates the result from Raghava n et al] 
(12010 1). who estimated that 5.9% of nearby dwarfs have com- 
panions with s > 2kAU. 

More critical to the analysis is whether the distribution of 
orbital separations is flat with log s or undergoes a signifi- 
cant decline. By limiting our analysis to subsets of pairs with 
larger proper motions, we can substantially reduce the optical 
contamination, although at the expense of sample size. Both 
the distributions of orbital separations s and the total compan- 
ion fraction derived from this restricted sample (bottom panel 
of Figure [TT]) are very similar to the results for the full sam- 
ple, and strongly support a decline with log s. We see only 
a slight reduction of the companion fraction with increasing 
^ cutoff: 3.6 ±0.3% for fi > lOOmas yr"' and 3.3 ±0.3% 
for /i > 150 mas yr"'. The slight reduction for stars with 
higher proper motions could be consistent with those stars 
being statistically older, and thus more susceptible to disrup- 
tion. In addition, the distribution in logp is also declining, 
as evidenced hy b = -0.29 in the model The slope is 
-0.29 X 1 .78 = -0.52 per decade (the range of p is 1 .78 dex), 
in excellent agreement with / = 1.5 found here by a different 
method. 

The distribution of physical companions for Solar-mass 
stars in the vicinity of the Sun in log s is thus declining faster 
than Opik's law: f(s) oc s~' with / ~ 1.5 (where / = 1 would 
correspond to Opik's law). This conclusion is valid at least 
over the 2kAU< s <20kAU range. A deficit of wide bina- 
ries at large separations (or periods) relative to the O pik's law 
was also d ocumented by Chaname & Gould] (l2005i) . and sug- 
gested by iLepine & Bong iorno (2007|) using a subset of our 
current data. This decline is supported by the current, ex- 
tended sample. Lepine & Bongiorno (2007) found a declin- 
ing distribution with / « 1.6 at 5 > 3 kAU; the slightly faster 
decline could however be the result of an additional selection 
criterion (A/i oc p"') which favors closer binaries. 

Figure [12] compares the empirical distribution of separa- 
tions with two possible analytical models: a simple power 
law, and the log-normal period distribution of Raghavan et alj 
C2010.) . For the latter, the median period 10^ d, logarithmic 
dispersion 2.28, and binary frequency 0.5 was converted to 
the distribution of semi-major axis by assuming total mass 
of 1.5M0 and total sample size A^ief = 5196. The agreement 
with the log-normal distribution is better than one might ex- 
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pect, given the fact that their sample does not overlap with the 
stars studied here and that the analysis methods are different. 

The frequency of 500-AU compa nions to FG dw arfs is 
higher, 13% per decade of separation (lTokovminll201 lb . This 
agrees with the distribution declining at long orbital peri- 
ods. Indeed, extrapolation from the first three bins (geo- 
metric mean s = 5.65 kAU, frequency 3.8% per decade) to 
s 500 AU with the s~^-^ power law predicts 3.4 times more 
companions, 13% per decade. In other words, the 5"'^ law 
seems to hold at separations less than 1 kAU. 

6. SUMMARY AND DISCUSSION 

We have searched for wide companions to Hipparcos stars 
in the SUPERBLINK proper motion catalog based on com- 
mon proper motion. The main results are: 

• The mass ratio of wide companions to nearby solar- 
mass dwarfs is found to be distributed uniformly to first 
order, with a possible deficit of companions with mass 
ratios q ~ 0.5. 

• The distribution of projected separation s declines as 
f(s) oc i"'^, faster than the Opik's law. 

• Total frequency of wide companions to solar-type 
dwarfs with 2<s< 64 kAU is 4.4 ±0.3%. The fre- 
quency in the range from 2 to 16kAU is 3.8 ±0.3% per 
decade of separation. Both numbers are lower limits 
because a fraction of wide companions is missed. 

In addition, our search has produced a catalog of 1413 CPM 
companions to Hipparcos stars. The present catalog gives the 



most complete census of wide CPM companions to nearby 
stars available to date. A decline in companion frequency at 
large s found here can be used to constrain gravitational per- 
turbations in the Galactic field. 

We found 21 potential WD companions to nearby stars (Ta- 
ble |2]l. Half of those should be true WDs, the rest are faint 
background stars. Indeed, 5 WD companions are already 
known in the literature, 6 look promising based on their sepa- 
ration and PM. 

Our cat alog contains new wide companions to exo-planet 
host stars (ISchneiderll201 lb . in addition to previously known 
CPM companions. Relevant data on two high-probability 
physical companions to HIP 17747 and HIP 90593 are listed 
in Table[3] the potential WD companion to HIP 6008 1 is listed 
in Table |2] More information on these candidates (spectral 
types, radial velocities) are necessary for their definitive attri- 
bution. 
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